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Abstract 

We explore the method of entanglement entropy applied to 2d black holes. We 
introduce a solvable model of a real scalar field with finite volume and lattice spacing 
in terms of N coupled mechanical oscillators and compute its entanglement entropy 
in many cases. The large N limit of this scheme, with finite lattice spacing, should 
reproduce the Bekenstein-Hawking formula for black hole entropy. 



1 Introduction 



Recently, remarkable progresses have been made in order to understand the microscopic 
origin of the thermodynamic behavior of black holes [l]-[2] ( for a review see [3] ). Let us 
mention the rigorous proof of black hole entropy achieved in superstring theory, using the 
D-brane technology [I]- [I]. Analogously in loop quantum gravity |8j black hole entropy has 
been identified with the logarithm of the total number of all different spin network states for 
a fixed eigenvalue of the area operator [9]. 

The derivations made in these candidates for quantum gravity strongly depend on the 
details of these theories. However it is generally believed that the proportionality of the 
black hole entropy with the horizon area is more universal and not dependent on the details 
of the theory. 

In this direction we can recall 't Hooft [lQj identification of black hole entropy with 
the statistical entropy of a thermal gas of quantum particles in the Hart le- Hawking state 
[TT]-[T2] with Dirichlet boundary conditions just outside the horizon ( the so called brick 
wall model ). He found, in the WKB approximation, apart from the expected volume 
dependent thermodynamic quantities describing hot fields in an almost flat space, additional 
surface contributions proportional to the horizon area. These extra contributions are also 
proportional to a~ 2 , where a is the proper altitude of the wall above the gravitational radius 
and they diverge in the a — > limit. These divergencies have been recognized as quantum 
corrections to the Bekenstein-Hawking formula that can be reabsorbed in the renormalization 
of the one loop effective gravitational lagrangian. For a specific choice of a, which is generally 
at the Planck scale, 't Hooft has been able to reproduce the Bekenstein-Hawking formula 
with the correct coefficient. 

The original formulation of the brick wall model has been corrected in [13]- [16]. Recently 
it has been reconsidered in the special case of (1 + 1) dimensions [17J, because it can be 
solved exactly without any approximation. 

A concurrent method has been developed in [H]-[20], i.e. the entanglement method, 
which we are going to reconsider in this article always in the solvable (1 + 1) dimensional 
case. This method is a strong candidate for a model independent derivation of the origin 
of black hole entropy, since it defines a statistical entropy measuring the loss of information 
induced by a spatial division of a system. 

This division distinguishes between the inaccessible quantum degrees of freedom made 
unphysical by the presence of black hole and the residual physical ones. It is expected that, 
independently from the details of the theory, entanglement entropy is proportional to the 
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area of the boundary of the spatial division of the system. 

The advantage of (1 + 1) dimensions is that there is a coordinate transformation making 
the background metric conformally flat, and therefore the quantization of a scalar field in 
the background of a black hole can be reduced to the free case. 

We point out that the natural spatial division which we adopt to compute black hole 
entropy is not separating the interior of black hole from the external region, but, since the 
coordinate transformation flattening the metric carries an extra length scale a, the wall 
proper altitude, is distinguishing between the interior of the brick wall and the external 
region to the wall. 

The quantization of a free scalar field is solved by introducing the proper modes, which 
are definitely non local with respect to the spatial division of the scalar field. Therefore 
the vacuum wave function carries correlations between the inaccessible degrees of freedom 
inside the brick wall and the residual physical ones. Such a correlation is the mechanism 
responsible for the production of a non trivial entropy, which is similar to the brick wall 
entropy calculated by 't Hooft. 

To simplify the calculation it is necessary quantizing the scalar field reduced to a chain 
of N coupled oscillators with periodic boundary conditions ( finite volume and finite lattice 
spacing ), i.e. a system with a finite number N of degrees of freedom. We are then able to 
carry out the entropy calculation in many cases of entanglement applying the general method 
of ref. [18]. In all these cases the entropy is finite and dependent only on the frequencies of 
the normal modes. 

However only in the N — > oo limit and finite lattice spacing it is possible to determine a 
surface entropy contribution ( proportional to the area of the horizon and diverging as a~ 2 , 
with a the lattice spacing, in four dimensions ). In this limit we have to solve exactly the 
spectrum of an integral operator, but we haven't found any better insight with respect to 
the qualitative features depicted in [18]- [19], so we leave its solution to a future study. 

2 The model 

The brick wall model [10] has shown that black hole entropy can be obtained in an heuris- 
tic way from the quantization of a scalar field on curved background extended to a finite 
temperature state, in which certain quantum degrees of freedom are lost because of the pres- 
ence of black hole. The identification of these degrees of freedom can be realized with the 
entanglement method that we are going to apply to the simplest case of (1 + 1) dimensions. 
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In this article we limit ourself to consider a scalar field immersed in a background metric 
{M, g} in 1 + 1 dimensions. The action 

S = j d 1+1 x^/\g~\ (^d^d v <f> - \m 2 <p 2 ^ (2.1) 
implies the Klein- Gordon equation on a curved space 

-^=d,(g^^\d,)<P + m 2 <p = (2.2) 



V\9\ 

From the action S (12. II) we can derive the corresponding energy-momentum tensor 

2 (99 1 

9, u (x) = = -<pd,d u <P + -g, u (g^dpd^ + m 2 <p 2 ) (2.3) 

y \g\ 9 \ x ) ^ 

We restrict the choice of the background metric to the case of a classical (1 + 1) dimen- 
sional black hole 

ds 2 = f(r)dt 2 - -i-rfr 2 (2.4) 

fin 

f(r) is an arbitrary function constrained to the following properties: 



lim f(r) = 1 

r— >oo 

/(r) has a simple zero at r = ro (2-5) 
where ro is the horizon with positive surface gravity 

ko = \f(r Q ) > (2.6) 
The classical equation of motion for the scalar field is 

' d 2 (t)-d r {f{r)d r <P\ +m 2 (f) 2 = (2.7) 



In the coordinates (t, y) with "tortoise" coordinate y defined by 



the black hole metric becomes conformally flat. Notice that /(ro) = and /'( r o) — 2&0 > 
implies that for r — > r at the leading order 

y ( r ) „ -L/n2fc (r - r ) /(r) ~ e 2k ™^ « 2£; (r - r ) (2.9) 
zfc 

and therefore, containing a logarithm, it introduces an auxiliary scale of length 1/2&Q. 
Moreover the coordinate transformation is singular near the horizon, hence in the y coor- 
dinates it doesn't make sense introducing a spatial division of the system coinciding with 
black hole horizon. 

Let r\ > ro be the radius of our alternative spatial division of the system, analogous to 
't Hooft brick wall, and let us calculate the proper altitude of the wall a over the horizon 
depending on Sr = r\ — r : 

f n i 1 
Jro 2 

In the coordinates (t, y) the Klein-Gordon equation (12.71) reduces to its flat form in the 
massless limit ( m — ) 

(d 2 t - d 2 )cj> = (2.11) 
and we have the problem of quantizing a free scalar field. 

This length scale a is also present in the brick wall model, where it can be proved to be 
of the Planck scale order and it allows to distinguish the regions 

y < y(r x ) r < r < n 

y>y(n) r > n (2.12) 

It is aim of our article to hypothesize that the degrees of freedom freezed by the presence of 
black hole are situated in the interior of the brick wall for y < y(ri), i.e. they are superficial, 
while the residual physical ones are always external to the brick wall for y > y{r\). In this 
sense we distinguish ourself from the literature that usually identifies them in the interior of 
black holes. 

We notice that in the y coordinates in two dimensions this spatial division between frozen 
and physical degrees of freedom seems the only possible choice. In greater dimensions or 
in the massive case it is not possible to avoid the presence of a background field and the 
analysis appears to be more intricate. 
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3 Scalar field as a finite chain 



Before quantizing the scalar field we must restrict ourself for simplicity to a model of it in 
which both the volume and the lattice spacing are finite. It is well known that the scalar 
field in (1 + 1) dimensions, defined on a lattice and with a finite volume, reduces to a system 
of iV coupled mechanical oscillators. 

Consider the hamiltonian 



1 N v 2 

H = -^[u\ + —{u A -u A+1 f + ulu A } (3.1) 



2^ L A • a 

A=l 

where the single oscillator variables u A ,A = 1,..,N are real and we impose periodic 
boundary conditions (N + 1 = 1). 

Such an hamiltonian can be diagonalized introducing the so-called normal coordinates 
through a discrete Fourier transform 



e ik(k)a A _ 2tt 

u a = / ^ — T^Tlk k = —k a A = a- A 
t! VN Na 

N e -ik(k)a A 

^ = E ~m Pk (3 - 2) 



or inversely 



^ e ~ik(k)a A 
^ e ik(k)a A 

» = g^T* 4 (3 ' 3) 

The inversion is obtained by recalling the discrete completeness relation 

N in(k—k')a 

£^^ = %' ( 3 - 4 ) 

n=l 

The reality condition for u A , A — 1, ...N is assured by the following necessary condition 
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p*k = p-k 



(3.5) 



By introducing the parameterization (qk,Pk) the Hamiltonian (13.11) can be recast in the 
following diagonal form 

k 

where ujk is given by the dispersion relation 

^ = -o 2 +4» 2 ^ (3.7) 

The information about the original coupling between the oscillators is now contained in 
the dispersion relation. The massless limit is reached by putting ujq = 0. 

Such a diagonal form of the Hamiltonian is ready for the quantization procedure of the 
coupled system of N oscillators. 

The canonical quantization rules for the variables ua 



[u a ,Ub] = i?i$AB 

[u A ,u B ] = [ua,u b ] = (3.8) 

can be extended to the normal coordinates qk, since they are obtained with their conjugate 
momenta from ua, Ua through an unitary transformation (13.31) 



[qk,Pk'] = ifrhk' 

[qk,qk>] = [pk,Pk>] = o (3.9) 

Now ua and ua are promoted to hermitian operators, u\ = ua, and the corresponding 
normal coordinates have the following hermicity properties 



Q-h 
P-k 



(3.10) 



The corresponding Hamiltonian operator, written in terms of the normal coordinates, 
takes the following form 

# ^E^t+^W) (3-ii) 



2 

k 

The quantization procedure is simplified if we introduce the creation and destruction 

t 

k 



operators at and a k 



dk = 



vm (uJkqk+tpl) 



4 = ^f— ^ k< il - iPk) (3-12) 



or inversely 



qk = ^k {ak + alk) 

Pk = i^={a- k -a\) (3.13) 
From the definition of a k , a\ we immediately get 

[a k ,av] = [al,a] t ,] = (3.14) 

Since a_ k ^ a\ we can state that there are 2N independent operators a k , a\. 
We can now express H in terms of the operators a k , a\ obtaining 

H = J2^k(ala k + ^) (3.15) 

k 

The state of lowest energy |0 > is determined by the condition 

a fc |0>=0 VA; (3.16) 

In the Schrodinger representation in which the q k are basic variables, this definition can 
be rewritten as a differential equation for the variables q k , solved by 
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/(?*) = foe 



(3.17) 



Such an eigenstate has eigenvalue E = J2k _ 2 iL - The g ener ic eigenvalues of are given 



by 




(3.18) 



k 

Introducing the normal modes we have passed from a description in terms of single 
oscillators ua to a description of collective coordinates qk- 

4 Entanglement entropy 

The quantization procedure is naturally related to the collective coordinates, while the spatial 
division of the system induced by the presence of the black hole is better described in terms 
of the single oscillators ua- Hence when we integrate on the single oscillator degrees of 
freedom contained inside the wall for y < y{r%), clarifying where is the loss of degrees of 
freedom induced by the black hole, there is a non trivial coupling given by the vacuum wave 
function between frozen and physical degrees of freedom. This correlation is the ultimate 
reason for the presence of a non trivial entropy. 

The vacuum wave function 




(4.1) 



rewritten in the single oscillator coordinates ua (13.21) becomes 




N 



A=l 





where 




uj k e 



ik(k)(a B 



«A) + <_> 5 ) 



(4.3) 



Imposing the normalization condition of the wave function we get 
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^o(ma) 



det- 



M 



7T 



e -\M AB u A u B 



(4.4) 



Mab is a positive definite quadratic form; this property will be crucial in the following. 
The density matrix for the vacuum is given by 



P = IV'o >< 

and in the Schrodinger representation it becomes 



(4.5) 



P{ 



u , u 



,B X 



det 



M 



71 



e 2 



\M AB (u A u B +u lA u' B ) 



(4.6) 



Now we take into account the presence of the black hole. We have said that, according 
to the hypothesis of section 2, in the region contained between the horizon and the brick 
wall, the quantum information is inaccessible. We can therefore obtain a reduced density 
matrix p rec i for the degrees of freedom outside the brick wall, integrating the single oscillator 
coordinates ua that are situated in the region y < y{r\): 



Pred (u a ,u' b )= I l[du a p(u a ,u a ,u' b ,u a ) 

J a 

_ 1 

u a u b +u' a u' b ) J Y^(l u a e i- M l3-yU' 3 u'l-M aa (u a +u' a )u '] ^ 7) 



where the matrix Mab is decomposed as 



M 



AB 



and the inverse matrix M AB [M AB Mbc — 8, 



M 



AB 





M ab 


M aa 




_ M ab 


M af3 


M BC = 5$) 






M ab 


M aa 




M ab 





(4- 



(4.9) 



If we introduce as in ref. [18] the following extra notations, i.e. M ab as inverse of M a & 
, M a b as inverse of M ab , M al3 as inverse of M a p e M a p as inverse of M a/3 , performing the 
quadratic integration we end up with the following formula 
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Pred(u a , v! ) 



det 



M, 



ab 



e 2 



\M ab (u a u b +u' a u' b ) ±M a ' 3 M aa M /3b (u+u') a (u+u , ) b 



71 



(4.10) 



where 



detM AB = detM ab detM c 



a 6 



(4.11) 



The integration is non trivial only because of the coupling given by M aa between the 
degrees of freedom inside the region y < y{r{) and those external to it (y > y(r{)). 

Using the identity 



M ab = M ab - M aa M a ^Mp h 



(4.12) 



and defining 



K ab = M ab H ah = M„ n .\r ! .\l 



rib 



(4.13) 



the resulting density matrix is 



pred{u a , u } 



det 



K. 



ab 



7T 



e 2 



\K ab (u a u b +u' a u' b ) 



e 4 



\H ab (u-u') a (u-u'f 



(4.14) 



In the special case of two degrees of freedom, of which one inside and the other outside 
the wall, the density matrix (I4.14p becomes 



p red (u,u') 



K 



71 



e 2 



K{x 2 +x' 2 )-\H{x-x') 2 



(4.15) 



The entropy of this mixed state depends only on the ratio 



(4.16) 



and on the accessory variable 
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„=l +5 -2 



1 I 1 
A ( A 



(4.17) 



through the formula 



S 



filnfi + (1 — n)ln(l — /i) 
1 — li 



(4.18) 



In the general case it is possible to extrapolate the general formula for the entropy by 
bringing the density matrix (14. 141) to the following diagonal form 



p red (u a ,u' b ) 



n 



Jl -\ e [-\(u a U a +u' a u' a )-\\a{u-u') a {u-u') a 



(4.19) 



where the A a are the eigenvalues of if ao , in the basis where K a b is the identity matrix, 
and in general they are recognizable as the eigenvalues of the composition of matrices 



A a b = (K-r c (H) cb 



(4.20) 



This is possible since at least one of K a b and H a b is associated to a positive definite 
quadratic form, with positive eigenvalues. Hence with a series of transformations K ab can 
be diagonalized simultaneously with the matrix H^. 



Finally, the entropy associated to a system of many degrees of freedom is given by 



Halnfia + (1 - fJ. a )ln(l - fj. a ) 



1 - Ma 



(4.21) 



where /i a is the positive solution of 



4/i a 



v 1 "Ma) 

and A a are the eigenvalues of the matrix A£\ 



(4.22) 



5 Entropy computation 

The natural application of the general theory outlined in the last section is the model of the 
scalar field reduced to iV degrees of freedom. We will be able to compute directly without 
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examining the reduced density matrix the entropy of such a system in terms of the frequencies 
ujk characterizing the normal modes of the field. 

In the simplest N = 2 case with a degree of freedom inside the brick wall and the other 
outside, we obtain from the general parameterization ( 14. 3D 



M 



AB 



— L0 X +UJ 2 UJ\+ <jJ 2 



T AB 1 ( UJ\ + UJ 2 U)i - U) 2 



M = —— - — - ~" (5.1; 

and therefore the matrices K and H reduce to pure numbers 

K= 2UlUJ2 H= ^- ^ (5.2) 
UJ1 + UJ2 2(0)1 + ^2) 

Finally the eigenvalue A ( I4.16P from which the mixed state entropy depends on is 



1 (UJI - Ld 2 



,2 



A = ^ (5.3) 

4 U>xU>2 

We can go on in generalizing this result to the case (1, N — 1), i.e. in which one degree 
of freedom is observable and N — 1 are the inaccessible ones. 

Let us take for example the case N = 3 with (1, 2), we obtain in particular 



Qk 
Qk 



SkAUA M AB = Re {S\ k to k S kB ) 



T ct 
u A a Ak 



(5.4) 



where Re implies taking the real part and we compute the transformation matrices S^a 



as 



(a a 2 1 \ 



V 1 1 1/ 

The specific form of Mab is 



ct L 
Ak ~ V3 



1 \ 



a a 
a a 

V 1 1 1/ 



.,2 x 



a 



_ 2-Ki 

e 3 



(5.5) 
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M AB 



' M ab 


M aa ' 


_ M ab 


M af3 



Pi 


P2 


P2 


P2 


Pi 


P2 


. pi 


P2 


Pi . 



where 



Pi 

P2 



-(wi +^2 + ^3) 



The inverse matrix satisfies the following properties 



M AB = 



M ab 






M a0 





$2 


$2 " 


S2 


Si 




. ^ 


S 2 


^ . 



where 



s = Pi + P2 

(Pi ~ P2){pi + 2p 2 ) 

5 2 = * 

{pi ~ P2){pi + 2p 2 ) 

In this case the matrix K is still a number 
while the matrix if is a number given by 

(P2 P2) ( Pi ~P2 \(P2 
p\~ p\ I ~P2 Pi I I P2 



H = M aa M a(S Mfs b 



and the eigenvalue A = ^ is given by 

A = 2pl = 4 (u 3 - ^ f 

(Pi ~ P2){Pl + 2p 2 ) 9 (^i+CJ 2 )^3 

It is also easy to study the case (2, 1): 
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M AB 



' M ab 


M aa ' 




M af3 



where p\ and p 2 are as in eq. (15.71) . 
Analogously the inverse matrix 



Pi 


P2 


P2 


P2 


Pi 


P2 


. 


P2 


Pi . 



(5.13) 



M 



AB 



M ab 




M ab 


M Q/3 



In this case the matrix K n h becomes 



while the matrix H a b is 



K, 



ah 



8 2 - 8 2 
°1 °2 



' 8 1 


S 2 


62 " 


82 


Si 


62 


. J 2 


62 


J l . 



5i -8 2 
-8 2 81 



(5.14) 



(5.15) 



H, 



ab 



P \l 1 1 

pi \ 1 1 

from which the matrix to be diagonalized is 



(5.16) 



H \h 61 ) \ 1 1 ) pi \ 1 1 

A£ is a matrix whose determinant is equal to zero while the trace 



Trk a h = 8 lPl -1 = 2- 



P5 



(pi -P2)(pi + 2p 2 ) 
from which its eigenvalues are recognizable as 



(5.17) 



(5.1* 



, 4(^3 

^1 — 7: 



tJl+tJg \ 2 
2 I 



A 9 



9 + u; 2 )^3 




(5.19) 



Therefore we can come to the conclusion that the cases (1,N — 1) and (N — 1, 1) are 
related by duality. At first sight in the case (N — 1,1) it seems that we have to diagonalize a 
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matrix of dimension (N—l), which is however singular since it contains N—2 null eigenvalues 
and only one that coincides with the case (1, N — 1). 

More generally we can guess the existence of a duality between the cases (k, N — k) 
and (N — k,k), with the first case containing k independent contributions to the total 
entropy, while the second one contains the same eigenvalues plus extra null contributions 
until completing N — k eigenvalues. 



6 N = 4 case 

Let us consider the case N — 4, that, apart from exhibiting the usual cases of the form 
(1, N — 1) or (N — 1, 1), introduces the new non trivial case (2, 2). 

Let us define the matrices Mab and M AB as 



M AB = Re{S\ k A km S rnB ) 
M AB M BC = 5 A 



(6.1) 



where the diagonal matrix A km has the following form 



A km 



UJ 2 



LJ 3 



\ 



(6.2) 



while the similitude matrix Ska is given by 



SkA ~ 2 



a 



a 



(a a 2 a 3 1 \ 



1 
1 

1 / 



^ 9 

a a a 



at 





/ 


a 3 


a 2 


a 


1 \ 


1 




a 2 


1 


a 2 


1 


2 




a 


a 2 


a 3 


1 




V 


1 


1 


1 


1 / 



Computing explicitly the matrix Mab we obtain 



S*S = 1 a 4 = 1 (6.3) 



M 



AB 



( Pi P2 P3 P2 \ 

P2 Pi P2 P3 

P3 P2 Pi P2 

V P2 P3 P2 Pi ) 



(6.4) 
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with the following definitions 



Pi = ^(cji + cj 2 + ^3 + ^4) 

P2 = ^(-^2 + ^4) 



while 



M AB = 



( S x S 2 S 3 S 2 \ 

5 2 Si S 2 S 3 

53 S 2 Si S 2 
\ 5 2 5 3 S 2 Si ) 



with the following definitions 



S 2 = - 
1 



Pi +P3 



1 



2 V(Pi+P3) 2 -4p| pi-p 3 

P2 

(pi + Ps) 2 - 4p| 



Pi +P3 



1 



2 V(pi+Ps) 2 -4p^ Pi-p 3/ 
Let us firstly discuss the (1, 3) case. 

The matrix M AB can be decomposed in terms of the following submatrices 



M 



AB 



M ab M aa 
M ab M afi 



while the matrix M is decomposed as 



M 



AB 



M ab M aa 



u 


P2 


P3 


P2 \ 


P2 


Pi 


P2 


P3 


P3 


P2 


Pi 


P2 


\ P2 


P3 


P2 


Pi J 


(Si 


S 2 


S3 


S2 \ 


S2 


Si 


S 2 


S3 


S3 


s 2 


Si 


S 2 


\ s 2 


S3 


S2 


Si) 
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Moreover M ab is the inverse of M ab = p\ and M ab is the inverse of M ab = 8±: 



M ab = 



M, 



Pi 



ah 



Si 



Finally M a/3 is the inverse of M a p 



M afB = — 



( Pi- Pi (P3 - Pl)p2 p\ - Plp3 ^ 

(P3 - Pl)p2 p\ - p\ (P3 - Pl)p2 
\ P\ ~ P1P3 (P3 - Pl)p2 pj ~ p\ ) 



with the determinant D defined by 



D = Pi{p\ - pl) - 2pl(p! - p 3 ) 



With these data we can compute the matrices 



K = M ab = ^ 



Si 



H = M U"M/„, = 



1 

D 



[p 2 3 (pi-p 2 3 ) + 2p 2 2 (pi + 2p 2 3 -3p lP3 )] 



Since D = S±(pi — p 3 ) 2 ((pi + p 3 ) 2 — ^p 2 ,) we can deduce that 



A 



" Pl) + IpWpI + 2P3 " 3pip 3 )] 

= (Pi + Ps)pl + 2(pi - 2p 3 )p\ 
~ (Pi-P3)((Pi+P 3 ) 2 -4p|) 



(6.10) 



(6.11) 



(6.12) 



(6.13) 



(6.14) 



Let us notice here that in general in order to compute the eigenvalues of the matrix A£ 
we can substitute M ac M c ^M^ a with — M aa and we can finally obtain 



a formula already encountered in the case (1, 2). 
Let us discuss the case (3, 1). 

In this case the decomposition of the matrix Mab is different: 



(6.15) 
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M 



AB 



M ab M aa 



and analogously for the matrix M AB 



M 



AB 



M ab M aa 

M ab M af3 



Pl 


P2 


P3 


P2 ^ 


P2 


Pi 


P2 


P3 


P3 


P2 


Pi 


P2 


\ P2 


P3 


P2 


Pi J 


/* 


6 2 


$3 




6 2 


Si 


$2 


S3 




52 


Si 


6 2 


\ 82 


S3 


62 


Si) 



The matrix M ab is the inverse of M ab 

( Pi - P\ (P3 - Pl)p2 p\ ~ PlP3 

(P3 ~ Pl)p2 pl ~ p\ (P3 ~ Pl)p2 

\ pi- pips {p3 - pi)p2 pi - p\ 



M ab = 



D(Pi) 



M ab is the inverse of M ab 



M ab = 



D(Si) 



( 51-51 (63-61)62 51 -5,5s \ 

(63 - 6^62 51 - 51 (63 - 5i)5 2 
V 5 2 2 - 5i53 {53-5i)5 2 51 - 51 ) 



M alS is the inverse of M a p = pi e M a p is the inverse of M alS = 5i. 

1 ~ 1 



M af3 



M, 



a0 



Pi 5 t 

From these we can recover that K ab = M ab is the inverse of M ab and therefore 



ah 



( 5i 5 2 5 3 \ 

5 2 5i 5 2 
V 63 6 2 5i J 



M an M^M pb = — 



( 



P2 



Pl 



P3 
\ P2 



(P2P3P2) 
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Being a projector it has two null eigenvalues and only one 7^ 0. 



A 



ab 



(K 
1 

Pi 



I 5 1 p 2 + 5 2 p 3 + 5 3 p 2 > 

25 2 p 2 + Sips 
\ 5ip 2 + 5 2 p 3 + 5 3 p 2 J 



(P2P3P2) 



(6.22) 



Let us notice that there is also the following simplification, due to eq. (16. 7j) 



{h + 5 3 )p 2 + 5 2 p3 = -5 2 pi 
25 2 p 2 + 5 1 p 3 = -5 3 pi 



(6.23) 



and therefore we obtain for the matrix A, 



ab 



( s 2 \ 



Kb = - S 3 (P2P3P2) 

Being of this form, all the matrix information is enclosed in its trace 



(6.24) 



TrA ab = -{2b 2 p 2 + d 3 p 3 ) = d 1 p 1 - 1 = — ; ■ — ^- (6.25) 

(Pi -Ps)((Pi +Pa) 2 -4p|) 

Being proportional to p\ and p\ the trace becomes null when all the frequencies are equal 
between them. K ab is a matrix with three eigenvalues, of which two null and one equal to 
the case (1,3). 

In general the cases (1, N — 1) and (N — 1, 1) are solved by the eigenvalue 



A = 8 x px - 1 



(6.26) 



where p\ = + uj 2 + .. + oj^) is the first entry of the matrix Mab and Si is the first 

entry of the inverse matrix M AB . 

Finally we discuss the case (2,2). 

In this case the matrix Mab is decomposed in the following submatrices 
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M 



AB 



M ab M aa 



Pi P2 
P2 Pi 



while the matrix M AB is decomposed as 



M 



AB 



M ab M aa 
M ab M afS 



P3 P2 
\ P2 p3 



( <$1 $2 

S 2 5i 



P3 P2 
P2 P3 



\ 



Pi P2 
P2 Pi ) 



(6.27) 



^3 S 2 
^2 ^3 



#3 2 
h 0~3 



5 2 \ 



51 5 2 

5 2 5i 



(6.28) 



/ 



The entropy calculation is in this case solved if we know the eigenvalues of the following 
matrix 



A a = _ M aa Mab = _ 




(6.29) 



The main reason why the eigenvalues of such a matrix A£ are real is because it can be 
rewritten as the product of two symmetric matrices of which at least one with all positive 
eigenvalues. This property assures that the two matrices can be simultaneously diagonalized. 

Let us notice that the matrix M ab can be diagonalized with the transformation 



P3 P2 
P2 P3 



P3 + P2 
P3 - P2 




(6.30) 



where a = b = from which 



M aa M ab 




S 3 + 5 2 
5 3 - 5 2 



P3 + P2 
p3 ~ P2 




(6.31) 



In this particular case the two matrices are simultaneously diagonalized by the same 
invertible matrix and we can directly compute the eigenvalues 



A± = -(fc ± S 2 )( P3 ± P2 ) = ^^4(P3 ± P2? 

\Pl ± P2> 

Therefore A± are real eigenvalues and the resulting entropy is real and positive. 



(6.32) 
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7 Continuum limit 



Until now we have limited ourself to the study of the entanglement of a scalar field reduced 
to a lattice of N point where the total entropy of the system is always finite. The continuum 
limit can be worked out by sending the volume to infinity and keeping the lattice spacing a 
finite, as in [IS]. In that case the entropy receives contributions from an integral operator 
admitting a spectrum of infinite eigenvalues accumulating around zero. The exact calculation 
goes beyond our possibilities and we are obliged to confirm the observations made in [18j-[19j, 
i.e. that it is possible to single out a finite surface entropy ( proportional to the boundary 
area ). In particular in |18j-[19j it is shown how to factor out in four dimensions an entropy 
per unit area: 

S/A = Ca~ 2 (7.1) 

where C is some constant and a is the momentum cutoff or lattice spacing. The pro- 
portionality constant C has the right order of magnitude if a is of the order of the Planck 
length l p . Although in general the system is not exactly solvable, it is possible to make some 
numerical estimate on the upper bound of this surface entropy. 

8 Discussion 

The entanglement method has been revealed effective to confirm the results obtained with 
the brick wall method. In this work we have firstly noticed a greater resemblance between 
the two methods clarifying that the inaccessible quantum degrees of freedom contributing 
to black hole entropy are situated inside the brick wall and not in the interior of black hole. 
We have pointed out that in two dimensions there is a coordinate transformation flattening 
the background metric and in the massless limit the quantization of the scalar field in the 
gravitational background can be reduced to the free field case. 

Then we have quantized the system reduced to a finite number N of degrees of freedom 
and computed in several examples their contribution to entanglement entropy finding some 
general result that may be useful in the large N limit. In these cases our model is an explicit 
application of the general method outlined in [18J. 

The ambitious goal would be to solve exactly the spectrum of the integral operator 
generalizing our results in the large N limit and to check if there is a contribution to entropy 
proportional to the area of the spatial division, but the technical difficulties are such that 
we are obliged to conclude our study here, postponing the large iV exact solution to a future 
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research. 
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